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Abstract 

W consider the problem of testing if a given matrix in the Hilbert space formula- 
tion of quantum mechanics or a function in the phase space formulation of quantum 
theory represent a quantum state. We propose several practical criteria to recognise 
states in these both versions of quantum physics. After minor modifications they 
can be applied to check positivity of any operators acting in a Hilbert space or 
positivity of any functions from an algebra with a Weyl type * - product. 

PACS numbers: 03.65.Ca, 03.65.Ta 



1 Introduction 



Mathematical structure of quantum mechanics has been formed from the two ingredients: 
observables and states. A notion of an observable is not unique and it depends on a 
chosen convention. In the Hilbert space formulation of quantum mechanics measured 
quantities are represented by linear operators. It is usually assumed that these operators 
are self -adjoint pQ, [2]. In the phase space formulation of quantum theory observables are 
represented by smooth real functions on a symplectic manifold [3] [6], but other functions 
can also be considered. 

On the contrary, a sense of a quantum state is precise. For example, in the Hilbert 
space formulation of quantum mechanics normalizable pure states are represented by 
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vectors of the unitary length from a Hilbert space H. Thus a set of all states is a convex 
set spanned by the pure states [7]. 

Formal construction of a space of states in quantum mechanics is based on a * - 
algebra A. The algebra A contains some subset of the set of observables. Let the symbol 
J 7 denote a set of linear functionals over the algebra A. Quantum states are represented 
by positive functionals / 6 J i.e. such that 

\/AeA f(A + A) > (1) 

and satisfying the normalization condition 

/(I) = 1- (2) 

By 1 we denote the unity of the algebra A. The expected value of an observable A G A 
in a state / G T equals 

(A) :=f(A). 

Detailed analysis of the algebraic attempt to quantum states can be seen in [8]. 

In the Hilbert space formulation of quantum theory the algebra A is the algebra 
B(H) of bounded operators defined on the whole Hilbert space "H. The functional action 
is determined by the formula 

f{A):=Tr{gA), A G B(H), 

where the symbol g denotes a density operator. More information about this topic can 
be found in [9]. 

The natural question arises, how to test whether a given functional represents a state. 
Definition 12.11 as well as geometric one are hardly applicable. Hence we propose several 
criteria for solving this problem. We consider the Hilbert space formulation and the phase 
space formulation of quantum mechanics. The most complex task is checking positivity. 
Our results referring to positivity are quite general and can be used e.g. in construction 
of a time operator [10J. 

The sign '*' plays a dual role in the text. It denotes both: an involution in the algebra 
A and a *- product. A square matrix C of a dimension dimH x dim'H is symbolised by 

1 dimW 

((fii\C\ipj) while an element of this matrix is represented as (<Pi\C\ipj) . 



2 Density operator 

Let us consider a quantum system on a separable Hilbert space %. Assume that {\<fij)}j^i 
is a complete set of ortho normal vectors in the Hilbert space H. We know only a probability 
of detecting the system in each of the states \ifj). To represent a state of the system we 
introduce a density operator. Following von Neumann [TTJ 
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Definition 2.1. The operator given by 



dim H 

q := u- lim V Pi |yj •) |, V j Pj > , V ^ = 1 

n— >dimTf ' » z — ' 

i=i i=i 



zs called a density operator. i£ac/i number Pj, j = 1,2, ... , dimH is the probability of 
observing the system in the state represented by a ket \<fij). If one of these numbers equals 
1 we say that the system is in a pure state. Otherwise the system is in a mixed state. 

The symbol u— denotes the uniform convergence of a sequence of operators. 

Definition 12. II presents a method of construction of a density operator g. Its equivalent 
formulation, in spirit of considerations presented in the Introduction, is the following. 

Definition 2.2. An operator g : H — >■ % is a density operator if it is: 



i. positive i.e. V|0) G H (<fi\g\4>) > 0, 
ii. self - adjoint g + = g and 
Hi. its trace Trp = 1. 



As we deal with operators acting in general in an infinite dimensional space, it seems 
to be necessary to comment on a notion of trace. For any bounded and positive operator 
A its trace is defined as [12] 

dim H 

Tri:= (fUlfi) 
i=i 

and is independent of the orthonormal basis chosen. So defined functional has properties 
analogous to the trace of a finite dimensional n x n matrix. Hence 

Tr(i + B) = Tri + TrS , V6 > Ti(bA) = b ■ TrA 

Moreover, the trace is invariant under any unitary transformation 

TriuAl/- 1 ) = Tri. 

From Defs. I2.ll andl 2.2l we can deduce several properties of the density operator. First 
of all we will see that the operator g is bounded. Indeed, by definition its domain is the 
whole Hilbert space H. Moreover, it is self - adjoint. Hence, from the Hellinger - Toeplitz 
theorem [12] we obtain that the density operator is bounded. 

As it can be seen from Def. 12.11 and from the spectral theorem, eigenvalues of the 
density operator are nonnegative and they do not exceed 1. Since the norm \\g\\ is the 
supremum of the eigenvalues of the operator g, it satisfies the inequality \\g\ \ < 1. 

From Def. I2.1l we easily calculate that Tig 2 < 1. Thus the density operator is a Hilbert 
- Schmidt operator. Its Hilbert - Schmidt norm, given by 

||£|| 2 := ^Tr(g+Q), 
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is not greater than 1. Moreover, as it is well known, ||£||2 = 1 if and only if the density 
operator represents a pure state. The Hilbert - Schmidt type operators constitute a * - 
ideal in the algebra B(H). 

Finally, as the eigenvalues of the density operator g are nonnegative, the density 
operator is a trace class operator and its trace norm 

| \g\ 1 1 := Tr \/ g* g — Tr\g\ = Tig = 1. 

The space of trace class operators B\{l-L) is also a * - ideal in the algebra B{%). For 
each trace class operator A the following estimation holds 

||i||<||i|| 2 <||i||i. (3) 

Since every density operator is a trace class operator, it is compact. For two arbitrary 
operators A e BiiTi) and B e B{%) the trace of the product is abelian 

Tr(AB) = Tr(BA). 

The density operator is positive. Hence for every operator A G B(9i) the mean value of 
the product obeys 

(AA + ) = Ti{gAA + ) > 0. 

In our paper we consider finite and infinite dimensional separable Hilbert spaces. In that 
second case two realisations: the set I 2 of one - column complex matrices and the space 
L 2 of square integrable functions over K are analysed. We discuss mixed and pure states. 



-r dim H 



Assume that a matrix representation (!fi\g\ipj) of an operator g is known. The 

kets {\(pi)}f™ H constitute an orthonormal basis of the Hilbert space %. We intend to 
settle whether this matrix represents a physical state of a quantum system. 



Let us consider a finite dimensional case dim'H < oo first. Then the matrix (<Pi\g\<Pj) 

completely determines the operator g, which is defined on the whole space H and is 
bounded. Therefore we use Def. 12.21 At the beginning we test if the matrix is symmetric 
Vy (<pi\g\(fj) = (<fj\g\(pi) and its trace ^fi™^ (<pi\g\ipi) = 1. Then calculating the prin- 

-i dim H 

cipal minors we check whether the matrix is positive. When the matrix (ipi\g\(fj) 
passes these three tests we conclude that it really represents a density operator. 



dim H 



For a pure state this procedure becomes simpler. A matrix ((pi\g\(fj) 
a pure state if it is symmetric, its trace equals 1 and the square of it is t 

dim "H 



dim H 



represents 

i 

re same matrix 



k=l 



The case of a separable but infinite dimensional Hilbert space is more elaborated. It 



may happen that in a given orthonormal basis the matrix (ifi\A\Lpj) 



of an 



J i 
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operator A exists but this matrix does not uniquely characterise this operator (compare 



Hence to know whether a matrix 



can represent a density operator, we 



(<Pi\Q\<Pj) 

propose to check first if the operator g is a Hilbert - Schmidt operator. As it was proved 
in p3], if 

oo 

^|(^|i|^-)| 2 <oo 

then the matrix [( | A | </}.,■)] represents an operator of the Hilbert - Schmidt type. Any 
Hilbert - Schmidt operator A is bounded and is defined on the whole space "H. Moreover, 
in this case the matrix (ipi\A\Lpj) completely characterises the operator A. In addition, 



if the matrix 



{ipi\A\ipj) is symmetric, we conclude that the operator A is self - adjoint. 



J 1 



Suppose that a matrix ((fii\g\ipj) is symmetric and of the Hilbert - Schmidt type. 



In the next step we have to test if the matrix 
We start from a few lemmas. 



(Vilfe-) 



represents a positive operator. 



Lemma 2.1. If a Hilbert - Schmidt operator A is positive, then in any matrix represen- 



tation 



(<Pi\A\(fj) all diagonal elements ((pi\A\<pi) are nonnegative. 



J 1 



Any Hilbert - Schmidt operator A is bounded. Every bounded positive operator is 
self - adjoint [12]. Moreover, each of its matrix representations determines this operator. 
Hence all the diagonal elements (ipi\A\fi) are real. From the definition of a positive 
operator for every % G M there is 



(<Pi\A\<Pi) > 0. 



Therefore we conclude that 



Corollary 2.1. A necessary condition for a Hilbert - Schmidt operator g to be a density 



operator is that in any matrix representation 
are nonnegative. 



(cpi\g\ipj) all diagonal elements {(pi\g\(pi) 



J 1 



Another observation is the following. 
Lemma 2.2. If g is a density operator then \ \g 



HI < 1. 



This statement is the straightforward consequence of the observation that the eigen- 
values of the operator g — 1 belong to the interval [—1,0]. As the norm of an arbitrary 
operator A is the supremum of the absolute value of its eigenvalues, the Lemma is proved. 



Lemma 2.3. If g is a density operator then the uniform limit of the sequence {(1 — g) n }^ =1 
of operators is the projective operator on the kernel of the density operator. 
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By Def. 12.11 eigenvalues of the operator (1 — g) are nonnegative and they do not 
exceed 1. Thus in the uniform limit we obtain a sum of 1-dimensional mutually orthogonal 
projective operators. This sum represents an orthogonal projection on the kernel of g. ■ 

The straightforward consequences of this lemma are two observations. 
Corollary 2.2. If g is a density operator then 

'u - lim (1 - g) n Y = u - lim (1 - g) n . 

\ n— s-oo / 71-+00 

Corollary 2.3. For every density operator g 



Tr ( (u — lim (1 

n— >oo 



e) n ) ■ g 



o. 



Therefore we propose the following criterion for a matrix (tpi\g\tpj) 
density operator g. 



J i 



to represent a 



Theorem 2.1. A matrix 

|2 



<V*|e|Vi) 



represents a quantum state iff: 



. V 1 < «, j < oo ((pi\g\tpj) = (VjlelVi); 



m. i/ie sequence 



(iPi\{i - Q) n \<pj) 



1 J n=l 



zs convergent in the norm, 



As it was shown, a density operator is a Hilbert - Schmidt operator and the trace of its 
square is not greater than 1. Moreover, it is self - adjoint. Lemma [2731 implies the property 
(in). By definition the trace of the density operator equals 1. 

The first condition says that the operator g is of a Hilbert - Schmidt type. As it is also 
symmetric, we see that it is self - adjoint. Thus from the spectral theorem it can be 
represented by a diagonal matrix and from (i) absolute values of its eigenvalues do not 



exceed 1. Therefore the matrix 
diagonal belong to the interval 

{<Pi\l - g\(pj) 



(<Pi\l — Q\<Pj) is also diagonal and the numbers on the 



J i 



0,2]. But if one or more elements of the diagonal matrix 



{<pt\0--$) n \<pj) 



n=l 



are greater than 1, the sequence of matrices 

is divergent. Hence the convergence of this sequence implies that all of eigenvalues of the 
operator g are nonnegative and the operator is positive. Thus its trace is well defined and 
if the condition (iv) holds, we see that g is really a density operator. ■ 

Another criterion is based on a notion of square root of an operator. 
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Lemma 2.4. [Wj Let A be a positive and self - adjoint linear operator defined on the 
Hilbert space H. Moreover, let \\A\\ < 1. Then there exists a unique positive and self - 
adjoint linear operator B such that B 2 = A. The operator B is a uniform limit of the 



series 



b = VA = Vi + (A- i) = i + \(A- i) - i(i - i) 2 + i(i - i) 3 + ... 

n\2 n v ' 



(4) 



A proof of this statement can be found in [12]. The series is convergent only for 
operators satisfying the estimation ||A|| < 1. Thus 



Theorem 2.2. A matrix 



((fi\g\(pj) represents a quantum state iff: 



M. V 1 < Z, j < OO (^1^1^-) = (^l^^i), 

ra. i/ie series 

- oo ^ - 



(¥>illlv>j) 



- oo ] 

, + 2 



(¥>tl£-ll¥>j) 



<¥>t|(£-i) 2 l¥>j> 



- oo ] 

i + 16 



<Vil(e-i) 3 lvj) 



+ 



-l) n+1 (2n-3)!! 
n!2™ 



+ . . . 



+ . . . 



(5) 



zs convergent in the norm to y (y>i|f?|<Pj) 

A proof of this theorem is analogous to the proof of Theorem 12. 11 The only difference 

r i oo 

refers to the series defining the square root of the matrix ((pi\g\ipj) . If the operator g 

is indeed a density operator, then it is positive and \ \g\ \ < 1. Thus its square root can be 
calculated as the sum of the series (jSJ). 

On the other hand, the series fl5]) defines the square root of an operator with the norm 
not greater than 1. Therefore if this series is convergent, the operator g is positive. 

By a slight modification of formula @ we obtain that the square root of a positive 
self - adjoint operator A satisfying ||A|| < 1 equals 



.. oo 



(-!)'+! (21 -3)!! 



1=2 



l\ 



2 l 



»-*)' ( I) Al - 

r=0 V 7 



(6) 



If the operator A is of a trace class so is each element of the series 
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Theorem 2.3. A matrix 



((fi\g\<Pj) represents a quantum state iff: 



*■ Eli K^lfe)! 2 < i. 



ii. V 1 < i,j < oo ((pi\g\<Pj) = (<Pj\Q\<Pi), 
Hi. the series 



1 r 



1 00 x — "V 

(Vi\Q 2 \(Pj)\ i +2^ 



(2/ -3)!! 



z-i 



2' 



D-D 



(7) 



«=2 r=0 

zs convergent to the operator g in the trace norm, 

The condition (m) says that £ = y^p. In fact it is sufficient to check this relationship 
in the sense of the uniform convergence of the series (J7J). However, this is a quite complex 
procedure. In this case it is much easier to test the trace convergence. 



Let us introduce a new symbol 



1 n 

7n ■= ^ + 



(-!)'+! (21 -3)!! 



Z=2 



/! 



E< 

r=0 



*2(I-r) 



When the sequence {(^ n — £)}J£Li is trace convergent then from the estimation (J3J) it is 
also uniform convergent. Thus the operator £ is positive. 

r 1 00 

On the other hand we can consider a diagonal form of the matrix (ipi\g\(pj) rep- 
resenting a density operator. For each index % the sequence {(<^i|£?n|y?i)}^i is conver- 
gent to (ipi\g\cpi). Moreover, each number ((pi\g n \(pi) is nonnegative and the sequence 
{( t Pi\Qn\^Pi)}'^ = i is growing for every i. Thus for every n the operator (g — g n ) is positive 
and its trace norm is simply the trace of the operator (g — g n ). The sequence of numbers 
{\\Qn — f?||i}5£=i tends to so the uniform convergence of the sequence {gn}^=i implies 
also its trace convergence to g. m 

The previous theorem leads to, perhaps, the most useful criterion formulated below. 



Theorem 2.4. A matrix 



((pi\g\(pj) represents a quantum state iff: 



*■ E«=il(^ilfe)| 2 < l > 



ii. V 1 < i,j < 00 ((pi\g\ipj) = ((pj\Q\(pi), 
Hi. the series 



1-1 



1=2 



r=0 
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w - E*=i (Vilely*) = !■ 



A density operator g is of a Hilbert - Schmidt type and self - adjoint. Moreover, it is 
positive, so A/tF = f?- As ||£|| < 1, the square root of the operator g is determined by the 



series 



Thus TT^J]P = 1. 



From the conditions (i) and (ii) we obtain that an operator g is of a Hilbert - Schmidt 
type and is self - adjoint. Its eigenvalues do not exceed 1. The square of a Hilbert - 
Schmidt operator is a trace class operator. Moreover, the set of trace class operators is an 
ideal in the space B(%). Thus the trace of each operator g 2n is well defined. The sum ([8]) 
equals Tvt,,/]P. If Try^ = Tr£, there must be yff? = g. Hence g is a positive operator. 
As its trace equals 1, it is a density operator. ■ 

The last two criteria presented in this section are based on a clever idea proposed by 
M. Wasiak [Tj 



Theorem 2.5. A matrix 



(¥>i|fe> 



represents a quantum state iff: 



ii. V 1 < z, j < oo ((pi\g\ipj) = ((pj\g\<fii), 
Hi. for every natural number n the sum 



E< 

fc=0 



)oo 
J2(<Pi\Q k+1 \Vi)>0, 
8=1 



The conditions (z), (ii) and (ii?) were considered before. Assuming that they are 

-i oo 

fulfilled, we can choose a basis {|<£i)}^i, in which the matrix (<pi\g\ipj) is diagonal. 
Nonnegative eigenvalues of the operator g are denoted by Xj, negative by — ^. 

Since K'aI^I'A)! 2 < 1) we see t na t & U < 1 and all z/i < 1. 



n /■ \ n / \ oo 

E(-i)*( £ W +I = B- 1 )' U ) E 

fc=0 ^ ' fc=0 ^ ' i=l 



a* • a* - Vi ■ (-l) k y k ) 



k=0 



-X 



i=l fc=0 



i=l i=l 

As n tends to oo, the first component of the sum goes to 0. On the contrary, the sum 
+ 2/i) n grows for n -)■ oo. ■ 

An obvious consequence of Theorem 12.51 and its proof is the following statement. 
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-i oo 



Theorem 2.6. A matrix {(pi\g\<Pj) represents a quantum state iff: 



*■ ES=iKvilelvj)l 




m- V 1 < i, j < oo (</?i|fe) = {<Pj\Q\<Pi), 



Hi. the limit 



fc=0 v 7 i=l 



iv - £i=l (<Pi\Q\<Pi) = L 

Presented criteria are of course valid not only for mixed but also for pure states. 
However, in this latter case the simplest test seems to be the following procedure. 



ii. V 1 < i, j < oo {(pi\g\(pj) = {<Pj\(}\(Pi), 

w - E*=i (<Pi\Q\<Pi) = L 

We begin with checking if the operator g is a Hilbert - Schmidt type, because only 
in this case the matrix (ipi\g\(pj) completely characterises the operator. Moreover, as we 
already mentioned, a symmetric Hilbert - Schmidt operator is self - adjoint. In the step 
(in) we test, whether the operator g is a projection. 

All of our previous considerations done for infinite dimensional matrices can be adapted 
for the Hilbert space L 2 of square integrable functions over R. Indeed, if functions 



Theorem 2.7. A matrix (<Pi\g\ipj) represents a pure quantum state iff: 



*■ £S=iKvilfe)l 2 = i 



f(x),g(x) e L 2 and g(y) = Af(x), 



we write 




Thus an operator A is symmetric if 



\fx,yeM A(x,y) = A(y,x). 



An operator A is a Hilbert- Schmidt operator if 
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A product of operators is calculated as 

(A-B)(x,z) = / A(x,y)B{y,z)dy 
Jr 

and a trace of an operator A can be found from the relation 

TrA = / A(x, x)dx. 
Jr 

3 Wigner function 

When the phase space formulation of quantum mechanics is considered, two fundamental 
elements: a phase space and a *- product must be taken into account. Examples of 
constructing a phase space can be found in [15] - [19] . In this section we restrict ourselves 
to problems, in which phase spaces are differentiable symplectic manifolds. 

There exists vast literature devoted to existence and construction of *- products. A 
systematic treatment can be found in [20J - [22J. On every symplectic manifold (M.,oo) 
there exists a nontrivial *- product [23] - [25]. An iterative construction of the natural 
*- product of the Weyl type has been proposed by Fedosov [26], [27J. We assume that a 
*- product is local and in its differential form is of the Weyl type. However, we see the 
necessity of applying an integral form of the *- product too. This postulate will become 
clear when we discuss some properties of Wigner eigenfunctions. 

An observable on a phase space (Ai,u) is any smooth real function on Ai being a 
formal series in the Dirac constant h 

oo 

C°°(M)M\ 3 A(q\ q 2n ) = KMq\ ■ ■ ■ , Q 2 1- (9) 

By q 1 , . . . , q 2n we denote local coordinates on the symplectic manifold (Ai,u). Each func- 
tion Ai(q l , . . . , q 2n ), dim M. = 2n, is a smooth real function on 

Following the scheme presented in the Introduction, we choose a *- algebra A. Our 
choice ensures that all smooth functions on (M, u), which are formal series in h and have 
compact supports, belong to A. We will return to the question of choice and structure of 
this algebra later on when a Stratonovich - Weyl correspondence will be discussed. 

The involution V is realised by the complex conjugation. A product in the algebra A 
is a Weyl type *- product. 

According to the general definition quoted in the Introduction, quantum states are 
positive linear functionals over the algebra A satisfying the normalization condition (T5]). 
As it was shown in [27], every such a functional / can be written in the following form 

VA(q\...,q 2n )eA 

f(A(q\...,q 2n ))= [ A(q 1 ,...,q 2n )*W(q 1 ,...,q 2n )t(q 1 ,...,q 2n )^ n = 
v 7 J M 
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/ W{q\..., q 2n )*A{ q \..., q 2n )t{ q \...,q 2n ) 



UJ 



■n 



(10) 



J M 



A real function 



oo 



C°°(M)[[h]}3t(q\...,q 2n ) 




i=0 



is called a trace density. An explicit construction of this function was proposed by 
Fedosov [28J. The trace density is determined by a symplectic curvature tensor and its 
derivatives. A function W^ 1 , . . . , q 2n ) contains information about the state. This function 
is called a Wigner function. The integral 



is often called a trace cLS IS cL classical counterpart of the trace of operator. 

A functional nowadays called a Wigner function on the phase space M 2n appeared in 
the literature eighty years ago [29]. Since in the case of a system with the phase space 
M 2n an explicit form of a mapping between operators and functions is known [3U], several 
properties of a Wigner function can be directly obtained from a density operator. They 
have been published in [3] [5], [3TJ and [32]. But even for systems modelled on the 
phase space M 2n the question of admissible states was considered only in (33] . Criteria of 
identification of pure states were proposed in [31] and [35] . 

In a general situation we do not know an explicit form of a mapping between operators 
acting in a Hilbert space and functions in a phase space. Yet, as both: the Hilbert 
space and the phase space formulation describe the same physical reality, there must be 
a correspondence between them. This correspondence implies an algebra isomorphism 
between the quantum algebra of operators An and the algebra of functions A on the 
phase space. This isomorphism is known as a Stratonovich - Weyl correspondence 
SW : Ay, —> A, [36]. Below we list properties of the Stratonovich - Weyl correspondence. 

1. The mapping is one - to - one. 

We stress a relationship between the algebras An and A. The choice of the algebra of 
operators An determines the choice of the algebra of functions A. On the other hand 
we know, that this choice is not unique and for different algebras An or equivalently 
A we obtain the same quantum mechanics. 

2. The SW correspondence is linear. Moreover, SW(A + ) = SW(A). 

3. SW(1) = 1 i.e. the constant function equal 1 on the whole symplectic manifold M. 
represents the identity operator. 

4. If an operator A is self - adjoint then SW{A) is a real function. 

5. The • - product of operators is represented by the *- multiplication of functions, 




(11) 



SW(A ■ B) = SW(A) * SW(B). 
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6. The trace of an operator equals 



Tri=-i— / SW(A)(q 1 , . . . , g 2n )t(g\ . . . , q 2n )oj n . (12) 



M 



The Stratonovich - Weyl correspondence establishes a correspondence between a density 
operator £ and a Wigner function W^g 1 , . . . , q 2n ) 

sw '((^ = wV ''''' q2n) ' (13) 



Applying the Stratonovich - Weyl mapping to a density operator we find properties of a 
Wigner function on an arbitrary symplectic manifold. 

Integration of a Wigner function yields 

W(q\...,q 2n )t(q\...,q 2n )u n = l. 

M 

This result is an immediate consequence of the fact that Tr^ = 1. Moreover, every Wigner 
function, as a counterpart of the self - adjoint operator, is real. 

Since a density operator is of a Hilbert - Schmidt type, the integral satisfies the 
estimation 

/; f! " 1 --' ! "»* ''>"^- 

The last equality holds only for pure states. The symbol W* m {q l , . . . ,q 2n ) denotes the 
ra-th power of the function W(q 1 , . . . , q 2n ) in the sense of the *- product. 

At the beginning of this paragraph we postulated (formula (JSj)) that observables are 
formal series in the Planck constant H. On the contrary, there exist Wigner functions 
containing arbitrarily great negative powers of the deformation parameter i.e. 

oo 

W(q\...,q 2n )= hW 3 {q\---A 2n ) (14) 

j=-oo 

and for arbitrary natural number m there is such a natural number s > m that 
W-si^q 1 , . . . , q 2n ) 7^ 0. Indeed, from the Stratonovich - Weyl correspondence and the fact 
that for a pure state its density operator is a projection operator, we see that for a pure 
state 

W* 2 (q\...,q 2 n = -^W( q \...,q 2 "). (15) 
Assume that W(q\ . . . , q 2n ) = Y,?=o MWjfa 1 , q 2n ). Thus 

oo 

W* 2 (q\...,q 2n )= £ h^ +k B k {WM\---A 2n ),WM\---A 2n ))- 

i,j,k=0 

The local bilinear operators £>&(-, •) define the *- product. 
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Thus from the relation ( IT5|) we see that Wo(q 1 , . . . , q 2n ) = 0. Using the mathematical 
induction we obtain that for every j there must be Wj(q 1 , . . . , q 2n ) = so we arrive to 
a contradiction. Therefore at least Wigner functions of pure states must be of the form 

(ED. 

Now we can show why in general an integral form of a *- product is required. If we 
calculate W* 2 (q 1 , . . . , q 2n ) as 

oo oo 
k=0 i,j—oo 

the term standing at h m , m G Z is an infinite sum 

£ B m ^(WM\ • • • , q 2n ), Wjiq 1 , q 2n )). 

i+j<m 

Unless the number of nonzero operators Bk(-,-) is finite, the sum above can be divergent. 
Thus in calculations involving Wigner functions we have to apply another, for example 
integral form of the *- product. Such a form is known for the Moyal product [6], [T8] . 

121]. 



The Stratonovich - Weyl correspondence between the Hilbert space formulation and 
the phase space formulation suggests that criteria of being a Wigner function should be 
based on a notion of trace. Indeed, in contrary to the uniform convergence of a sequence of 
operators, which does not have an easy applicable counterpart in the phase space version 
of quantum mechanics, the trace is simply represented by the integral ( II 2p . Therefore we 
can transform Theorems I2.4| 12.51 and 12.61 into their phase space versions. 

We start from a counterpart of Theorem 12.41 

Theorem 3.8. A function W(q 1 , . . . ,q 2n ) defined on a symplectic manifold (A4,u) is a 
Wigner function iff: 

1. j M W* 2 (q\ . . . , q 2n ) t(q\ . . . , q 2n )u n < 

2. the function is real, 

3. the series 

(-l)'+i(2Z-3)!! 



x 



( 1 ) (2nh) 2 ^-V [ W * 2 ^\q\ . . . , q 2n ) t(q\ . . . , g 2 >™ = 7^ 



i-i 

x 

r=0 



(16) 



I $ M W{q\...,q 2n )t(q\...,q 2n )uj n 
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An analog of Theorem 12.51 can be formulated in the form 

Theorem 3.9. A function W{q Y , . . . ,q 2n ) defined on a symplectic manifold (j\4,u) is a 
Wigner function of a quantum state iff: 

1. j M W* 2 (q\ . . . , q 2n ) t{q\ . . . , q 2n )u n < jafer, 

2. the function is real, 

3. for every natural number m the sum 

( ™ \ (27th) nk [ W< k+1 \q\ q 2n ) t{q\ q 2n )u n > 0, (17) 
fc£ V k ' J M 



I j M W{q\...,q 2n )t{q\...,q 2n )u n = l. 

Finally, a modified Theorem 12.61 reads as follows 

Theorem 3.10. A function W(q 1 , . . . , q 2n ) defined on a symplectic manifold (j\4,u) rep- 
resents a quantum state iff: 

1. j M W* 2 (q\ . . . , q 2n ) t(q\ . . . , g 2 >« < 

2. the function is real, 

3. the limit 

lim V(-l) fc f ™}(2irh) nk [ W< k+1 \q\...,q 2n )t{q\...,q 2n )u n = 0, (18) 
m ^°° t^o \ k J JM 



I J M W(q\...,q 2n )t(q\... J q 2n )u J n = l. 

Proofs of the presented criteria follow directly from their Hilbert space counterparts 
and the Stratonovich - Weyl correspondence. 

An identification method of pure states is based on the following statement. 

Theorem 3.11. A function W{q 1 , . . . , q 2n ) defined on a symplectic manifold (Ai,u) rep- 
resents a pure quantum state iff: 

2n\, ,n 1 

UJ 



1. f M W* 2 (q\...,q 2 nt(q\...,q - ^ 

2. the function is real, 

3. W* 2 (q\ . . . ,q 2n ) = j^Wiq 1 , . . . ,q 2n ) 

I j M W(q\...,q 2n )t(q\...,q 2n )u n = l. 
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We emphasise that the main problem with testing Wigner functions is calculating their 
*- product. As we said before, a differential form of this product can be useless in this 
case. 

Example 

As an illustration of the presented criteria we examine a function considered by 
Tatarskij in his distinguished paper [3]. The function 

W(p, q) = ^W (p, q) + ^Wxip, q) - ^W 2 (p, q) (19) 

is defined on the phase space R 2 . By Wo(p, q), Wi(p, q) and W 2 (p, q) we denote the Wigner 
functions of mutually orthogonal states i.e. 

/ Wi(p,q) *Wj(p, q)dpdq = — — 5^ . 

The function W(p, q) is not a Wigner function, because one of its eigenvalues is negative. 
Let us retrieve this observation by application of our three tests. 

As it can be easily calculated 

W* m (p, q) = ^^ZT ( (|) m Woip, q) + (f) m W 1 (p, q) + " W 2 (p, q)j . 

The function W(p, q) is real, the integral J R2 W(p, q)dpdq = 1. Moreover, 

J R2 W* 2 (p,q)dpdq = ^ so the function W(p, q) satisfies the conditions (i), (ii) and (iv) 

of Theorems ESI EH and IXTU1 

However, for the function f lT9|) the sum f lT6|) equals | • i so from Theorem 13.81 this 
function is not a Wigner function. 

Now in Theorem 13.91 we see that the sum (fT7|) equals 2nh for m = 0, for m = 1 and 
— | ■ 2nh for m = 2. Therefore after taking three steps we conclude that W(p, q) does not 
represent any state. 

Finally, the limit f|T8|) in Criterion 13.101 equals — oo so the tested function obviously 
cannot be a Wigner function. 

After a slight modification this example can be considered on any 2-D symplectic 
manifold. 
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